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We study the manifestly exotic tetraquark Di (cwsirf) and the scalar tetraquark fo(udud) in 
SU(3) C anisotropic quenched lattice QCD. We adopt the (9(a)-improved Wilson (clover) fermion 
at j3 = 5.75 on 12 3 x 96 with renormalized anisotropy a s /a t = 4, and investigate the correlator of 
the four-quark (4Q) system, cusd, with various quark masses including the idealized SU(4) / case. 
For fo(udud) etc., we only consider connected diagrams at the quenched level, i.e., the tetraquark 
fo(udud) is identical with D^ + (cusd) in the idealized SU(4)^ case. First, for comparison, we 
study the lowest qq scalar meson, and find that it has a large mass of about 1.37GeV after chiral 
extrapolation, which corresponds to /q(1370). Second, we investigate the lowest 4Q state in the 
spatial periodic boundary condition, and find that it is just a scattering state of two pseudoscalar 
mesons, as is expected. Third, to extract spatially-localized 4Q resonances, we use the Hybrid 
Boundary Condition (HBC) method, where anti-periodic and periodic boundary conditions are 
imposed on quarks (c,u) and antiquarks (s,d), respectively. By applying the HBC on a finite- 
volume lattice, the threshold of the two-meson scattering state is raised up, while the mass of a 
compact 4Q resonance is almost unchanged. In HBC, we find a nontrivial 4Q resonance state 
about 100 MeV below the two-meson threshold in some quark-mass region. Its chiral behavior 
largely differs from a two-meson scattering state. The scalar tetraquark fo(udud) is found to have 
the mass of about l.lGeV after chiral extrapolation, and seems to correspond to /o(980). Then, 
the manifestly exotic tetraquark D^(cusd) would exist around lGeV in the idealized SU(4)/ 
chiral limit. Finally, MEM analysis is applied to obtain the spectral function of the 4Q system. 
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1. Introduction 

The tetraquark [1-7] is an interesting subject in hadron physics. The tetraquark picture was first 
proposed for scalar mesons [[|]. Recently, charmed tetraquark candidates were discovered [|I], Qj. 

1.1 Experimental discoveries of tetraquark candidates 

In recent years, various candidates of multi-quark hadrons have been experimentally observed 
[g, and theoretically studied [4-9]. For instance, D+(2317) [0], D+(2460), X(3872) [@| and 
Y(4260) are expected to be tetraquark states [Q, ^] from the consideration of their mass, narrow 
decay width and decay mode. As the unusual features of X(3872), its mass is rather close to the 
threshold of D°(cw) and D°*(mc), and its decay width is very narrow as T < 2.3MeV (90 % C.L.). 
These facts seem to indicate that X(3872) is a tetraquark ^ or a molecular state of D°(cii) and 
D°*(mc) rather than an excited state of a cc system [g]. Also, D^ } (2317) is conjectured to be not 
a simple meson of cs but a tetraquark from its small mass and its narrow decay width. In fact, 
simple quark-model assignment for D^ } (2317) is cs( 3 Po), but almost all theoretical predictions for 
P-wave cs states fail to reproduce its small mass [§, [l(J. Furthermore, there is a puzzle in the decay 
mode of 0^(2317): CLEO collaboration ^ shows that the dominant decay mode is Dfn {) and its 
contribution is much larger than the radiative decay as 



r(D+^D*+y) 
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^< 0.059. (1.1) 



If D^(2317) is isoscalar (1=0) such as cs, this decay pattern implies that the main decay mode 
is isospin breaking process of 0(e 4 ) in QED, which is rather anomalous. To explain this decay 
pattern, Terasaki proposed 1=1 possibility of D^(2317) [0], which means the tetraquark picture 
for D^(2317) and leads to as the isospin partner. 

Then, we examine manifestly exotic tetraquark D+ + (cusd) with C = +l,5 = +l and I$ = +l 
in anisotropic quenched lattice QCD. In this paper, we mainly report the idealized SU(4)y case. 

1.2 Tetraquarks in the light-quark sector 

There are five ++ isoscalar mesons below 2GeV: / (400-1200), / (980), /o(1370), /o(1500) 
and /o(1710). Among them, /o(1370) is considered as the lowest qq scalar meson in the quark 
model JTTj] . For, in the quark model, the lowest qq scalar meson is 3 Pq, and therefore it turns to be 
rather heavy. [/o(1500) and /o(1710) are expected to be the lowest scalar glueball or an ss scalar 
meson.] Then, what are the two light scalar mesons, /o(400-1200) and /o(980)? This is the "scalar 
meson puzzle", which is unsolved even at present. As the possible answer, in 1977, Jaffe proposed 
tetraquark (qqqq) assignment for low-lying scalar mesons such as / (980) and a (980) [||. 

Then, we examine the scalar tetraquark fo(udud) in the light-quark sector in anisotropic 
quenched lattice QCD. In this paper, we only consider connected diagrams at the quenched level. 
Note here that, if only the connected diagram is taken account at the quenched level, the tetraquark 
fo {udud) is identical with D^, + (cusd) in the idealized SU(4) f case. 

2. Lattice QCD 

We investigate the manifestly exotic tetraquark D^ + (cusd) and the scalar tetraquark f^udud) 
in anisotropic quenched lattice QCD mainly in the idealized SU(4) * case. 
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2.1 Diquark-antidiquark-type interpolating field operator 

For the tetraquark D^, + (cwsJ) with spin 7 = and isospin 1 = 1, we adopt the diquark- 
antidiquark-type (non-two-meson-type) interpolating field, 

O = e abc e ade {c T b Cy 5 u c ){s d Cy 5 d^), (2.1) 

to extract the 4Q resonance rather than two-meson scattering states. Here, the roman indices de- 
notes color indices, and C = 7472 the charge conjugation matrix. In this diquark-antidiquark-type 
operator, the overlap with two-meson scattering states is expected to be small. The energy of the 
low-lying state is extracted from the temporal 4Q correlator G(t) = hY,x(O(t,x)O^(0,0)), where 
the total momentum of the 4Q system is projected to be zero. To reduce highly-excited-state com- 
ponents, we use a spatially-extended source of Gaussian-type with p=0.4fm in the Coulomb gauge 



12]. We impose the Dirichlet boundary condition in the temporal direction [Bp. We analyze the 



effective mass m e ff(?) obtained from the 4Q correlator G(t) as m e ff(/) = \n{G(t) / G(t + 1)}. 

2.2 Anisotropic lattice QCD 

To get detailed information on the temporal behavior of the 4Q correlator G(t), we adopt 
anisotropic lattice QCD with the standard plaquette action and the 0(a , )-improved Wilson (clover) 



fermion at j8 = 2N c /g 2 = 5.75 on 12 3 x 96 with renormalized anisotropy a s /a t = 4 []|, [l2|]. This 
is anisotropic version of the Fermilab action [13]. The scale is set by the Sommer scale 1 = 



395MeV, which leads to the spatial/temporal lattice spacing as aj ~ l.lOGeV (i.e., a s ~ 0.18fm) 
and aj x ~ 4.40GeV (i.e., a t ~ 0.045 fm). The spatial lattice size is L = \2a s ~ 2.15fm. We use 1827 
gauge configurations, which are picked up every 500 sweeps after the thermalization of 10,000 
sweeps. On the happing parameter, we take fc=0. 1210-0. 1240 for light quarks and ff=0. 1120 for 
the real charm-quark mass. We summarize the lattice parameters and related quantities in Table 1 . 

Table 1: The lattice parameters and related quantities. They are almost the same as in Ref.[K|. 



p 


lattice size 


a s 1 a, 1 7g 


u s u t 


K 


5.75 


12 3 x 96 


l.lOGeV 4.40GeV 3.2552 


0.7620 0.9871 


0.1210(10)0.1240, 0.1120 



2.3 Hybrid Boundary Condition (HBC) method 

To extract low-lying 4Q resonances, we use the "hybrid boundary condition" (HBC) []8|] where 
we impose the the anti-periodic boundary condition for quarks (c, u) and the periodic boundary 
condition for antiquarks (s,d), as shown in Table 2. By applying the HBC on a finite lattice with 
L 3 , the two-meson threshold is raised up, while the mass of a compact 4Q resonance is almost 
unchanged. Then, the 4Q resonance may become visible as a low-lying state in HBC, if it exists. 



Table 2: The hybrid boundary condition (HBC) to raise up the threshold of two-meson scattering states. 





c, u 


s, d 


gg-meson 


2-meson threshold 


tetraquark (qqqq) 


PBC 


periodic 


periodic 


periodic 


mi +W2 


periodic 


HBC 


anti-periodic 


periodic 


anti-periodic 


Lt=i,2 v^+^pL 


periodic 
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Mqq [GeV] qq scalar meson 




mi [GeV 2 ] 

Figure 1: The lowest mass Mqq of the qq scalar 
meson plotted against m 2 n . The qq scalar meson 
mass Mqq is extracted from the connected diagram 
in quenched lattice QCD. 



ITleff(t) 




PBC 

HBG »■■*- 
threshold - - - ■ 



t 

Figure 2: The effective mass plot of 4Q system in 
PBC and HBC. The solid lines denote the plateau 
of m e fi(t), and the the dashed lines the two-meson 
threshold. The Dirichlet boundary condition is 
imposed in the temporal direction. 



For a compact 4Q quasi-bound state of cusd, since it contains even number of quarks, it obeys 
periodic boundary condition (PBC) in HBC, and therefore its energy in HBC is almost the same as 
that in PBC For a two-meson scattering state, both mesons have non-zero relative momentum 
Pmm = (if jifjif)' i- e -> \Pmm\ = V^Tt/L in HBC, while they can take zero relative momentum 
Pmin = in PBC. In fact, the two-meson threshold m\ +ni2 (mi, ni2'. meson masses) in PBC is 
raised up in HBC as ~ Y,k=i 2 \l m \ + Pmm w i tn \Pmm\ = V^ti/L ~ 0.5GeV for L ~ 2.15fm. 



3. Lattice QCD results for tetraquarks 

First, for comparison, we calculate the qq scalar meson mass in anisotropic quenched lattice 
QCD. Here, we only consider connected diagrams. We show in Fig.l the lowest qq scalar meson 
mass plotted against m\. Our quenched lattice QCD indicates that the lightest qq scalar meson has 
a large mass about 1.37 GeV after the chiral extrapolation. Thus, the lightest qq scalar meson cor- 
responds to /o(1370) or ao(1450), which is consistent with the quark-model assignment for scalar 
mesons. [Note that, if disconnected diagrams are dropped off, isoscalar qq mesons degenerate 
isovector qq mesons in quenched QCD, e.g., m p = m w and m(/n) = m(ao).] 
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Lattice data (PBC) ' — • — i 






two-meson threshold (PBC) 
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Figure 3: The lowest energy M(4Q) of the 4Q 
system in PBC plotted against irv\. The dotted 
curve denotes the two-meson threshold of 2m n . 




0.2 0.4 0.6 0.8 1 

mi [GeV 2 ] 



Figure 4: The lowest energy M(4Q) of the 4Q 
system in HBC v.s. irv\. The dotted curve is the 
two-meson threshold of ym| + 3n 2 /L 2 in HBC. 
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Second, we investigate the 4Q system of cusd with various quark masses in PBC and in HBC. 
Here, we show the idealized SU(4)/ case, where D^ + (cusd) is identical with fo(udud) as far as 
only connected diagrams is considered. Figure 2 shows the effective mass m e ff(?) of cusd (or udud) 
in PBC and in HBC for fc=0.1240. In PBC, the plateau of m e s(t) agrees with the two-pseudoscalar- 
meson threshold, i.e., 2m n . As shown in Fig. 3, the lowest energy of the 4Q system in PBC plotted 
is just a %-% scattering state, as is expected. [Note that this lowest 4Q level is raised up in HBC, 
while a compact 4Q resonance does is almost unchanged in HBC] 

Third, we analyze the lowest 4Q state in HBC. In HBC, there appears a plateau below the n-Tl 
threshold by about lOOMeV as shown in Fig.3. Figure 4 shows the lowest energy of the 4Q system 
in HBC plotted against mr % . Note that its chiral behavior is largely different from the %-% scattering 
state. Thus, this 4Q state is considered to be a nontrivial 4Q resonance instead of a two-meson 
scattering state. After the simple linear chiral extrapolation, the mass of the 4Q state is estimated to 
be about l.lGeV, which seems to correspond to /o(980) or ao(980) in the light-quark scalar sector. 
This also indicates that the manifestly exotic tetraquark D+ + (cusd) would exist around lGeV in 
the idealized SU(4) f chiral limit. 

4. MEM analysis for point-point 4Q correlation 



Finally, using the maximum entropy method (MEM) [|14J, we extract the spectral function 
A (ft)) from the temporal correlator G(t), which satisfies 

,-(j3-T)£0 



G(T) 



dcoK(r,co)A(co), K(r,(o) = 



e~™ + e- 



l- e -Pco 



P =N t a t . 



(4.1) 



From the lattice data of the 4Q correlator G(t), we obtain the spectral function A(co) of the 4Q sys- 
tem with MEM. At present, we perform the MEM analysis for point-source point-sink 4Q correlator 
G(t) both in PBC and in HBC, although this correlator includes large excited-state contamination. 
For the point-point correlator G(t), the lowest-order default function m{co) is calculated as 

4N r 



m{(0) 



2 8 r(5)r(6)?r 6 



CO" 



(4.2) 



Unlike the previous section, we here use k = 0.1240 for u, d, s-quarks, and K = 0.1 120 for c-quark, 
and impose the periodic boundary condition in temporal direction. Figures 5 and 6 show the MEM 
analysis for the 4Q correlator G(t) and the obtained spectral function A (ft)) for the 4Q system. 



I PBC : Correlator 
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I PBC : Spectral hunction 1 
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a I GeV ] 



Figure 5: The 4Q corelator of point-source point- 
sink correlator G(t) in PBC. The symbols denote Figure 6: The spectral function A(a>) obtained 
the lattice data. The solid curve is the temporal with MEM from the point-source point-sink 4Q 
correlator reconstructed by A(co) in Fig. 6. correlator in PBC. 
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The lattice data of 4Q correlator G(t) can be well reconstructed with the spectral function A (ft)) ob- 
tained by MEM through Eq.(|4.l[). [In this point-point correlator G(t), highly-excited contamination 
is so large that low-lying structure of A((o) cannot be seen clearly. To clarify the low-lying structure 
of the 4Q spectrum, we are performing the MEM analysis for extended-source extended-sink 4Q 
correlators, where excited-state components are reduced.] 

5. Summary and conclusions 

We have studied the manifestly exotic tetraquark D^(cusd) and the scalar tetraquark fo(udud) 
in SU(3) C anisotropic quenched lattice QCD mainly in the idealized SU(4)y case. First, we have 
found that the lowest qq scalar meson has a large mass of about 1.37GeV after chiral extrapolation, 
which corresponds to /o(1370). Second, we have found that the lowest 4Q state is just a scattering 
state of two pseudoscalar mesons. Third, using the Hybrid Boundary Condition (HBC) method, 
we have found a nontrivial 4Q resonance state whose chiral behavior is different from a two-meson 
scattering state. The scalar tetraquark f^udud) is found to have the mass of about 1.1 GeV after 
chiral extrapolation, and seems to correspond to /b(980). Then, the manifestly exotic tetraquark 
D+ Q + (cusd) is expected to exist around lGeV in the idealized SU(4)y chiral limit. 
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